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摘要 : 

[目的 ] 通过 生物 信息 学 方法 探讨 METTL16 在 泛 癌 中 的 诊断 价值 及 肿瘤 预后 评估 的 作用 。 

[方法 ] 从 多 个 相关 数据 库 中 获取 各 个 肿瘤 及 正常 组 织 中 的 基因 数据 , 分 析 了 METTL16 在 
不 同 癌 症 中 的 表达 的 表达 水 平 ， 其 中 包括 癌症 细胞 系 百科 全 书 (CCLE )、 癌 症 基 
因 组 图 谱 CTCGAO 和 基因 型 -组 织 表 达 项 目 (GTEx )。 采 用 Cox 回归 分 析 和 
Kaplan-Meier 分 析 来 评估 METTL16 表达 水 平 高 低 与 肿瘤 预后 之 间 的 关系 。 使 用 肿 
瘤 免 疫 评估 资源 (TIMER) 来 评估 METTL16 与 癌症 中 的 免疫 细胞 浸润 的 相关 性 。 
同时 ， 还 分 别 分 析 了 METTL16 的 表达 水 平 与 常见 的 免疫 调节 因子 、DNA 甲 基 转 
移 酶 (DNMTs)、\ 错 配 修复 CMMRs)、 微 卫星 不 稳定 性 (MSI) 和 肿瘤 突变 负荷 (TMB) 
之 间 的 关系 。 此 外 ， 通 过 基因 和 集 富 集 分 析 (GSEA) 来 识别 癌症 中 的 关键 通路 。 

[结果 ] METTL16 在 多 种 癌症 中 异常 表达 ， 生 存 分 析 表 明 其 表达 水 平 与 许多 癌症 患者 的 临 
床 预 后 相关 。METTL16 的 表达 水 平 与 结肠 癌 (COAD)、 肾 透明 细胞 癌 (KIRC)、 
肝癌 (LIHC)〉 中 的 的 六 种 免疫 温润 细胞 (B 细胞 、CD4+T 细 胞 、CD8+T 细 胞 、 
树 突 状 细胞 、 巨 哈 细 胞 和 中 性 粒 细胞 ) 的 表达 水 平 呈正 相关 。 此 外 ，METTL16 表 
达 与 多 种 癌症 中 的 免疫 检查 点 基因 表达 、MMR、DNA 甲 基 转 移 酶 水 平 、TMB 和 
MSI 显 着 相关 。 同 时 ，GSEA 还 揭示 了 METTL16 参与 了 肿瘤 和 免疫 相关 信和 号 通 
路 的 调控 。 

[结论 ] METTL16 与 肿瘤 预后 和 肿瘤 免疫 浸润 密切 相关 ， 并 在 多 种 肿瘤 的 发 生发 展 中 起 着 
重要 作用 。METTL16 可 能 成 为 一 种 很 有 前 途 的 各 种 肿瘤 预后 生物 标志 物 ， 并 为 天 
来 的 肿瘤 免疫 研究 和 治疗 策略 提供 新 的 思路 。 
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[Objective] To explore the diagnostic value and prognosis of METTLI6 (METTLI16) in 
pan-carcinoma by bioinformatics method. 


[Methods| The gene data of each tumor and normal tissue were obtained from several relevant 


databases, including the Cancer Cell Line Encyclopedia (CCLE), the Cancer Genome Atlas 
(TCGA) and the Genotype-Tissue Expression Project (GTEx). The expression levels of METTL16 
in different cancers were analyzed.Cox regression analysis and Kaplan-Meier analysis were used 
to evaluate the relationship between METTLI6 expression level and tumor prognosis.The Tumor 
Immune Evaluation Resource (TIMER) was used to assess the association of METTLI6 with 
immune cell infiltration in cancer. Meanwhile, the relationship between the expression level of 
METTL16 and common immunoregulatory factors, DNA methyltransferase (DNMTs), mismatch 
repair (MMRs), microsatellite instability (MSI) and tumor mutation load (TMB) was also 
analyzed.In addition, gene set enrichment analysis (GSEA) was used to identify key pathways in 
cancer. 

[Results] METTLI16 was abnormally expressed in a variety of cancers, and survival analysis 
showed that its expression was associated with clinical outcomes in many cancer patients.The 
expression level of METTL16 was positively correlated with the expression levels of six types of 
immunoinfiltrating cells (B cells, CD4 + T cells, CD8 + T cells, dendritic cells, macrophages and 
neutrophils) in colon cancer (COAD), kidney clear cell carcinoma (KIRC) and liver cancer 
(LIHC).In addition, METTL16 expression was significantly associated with immune checkpoint 
gene expression, MMR, DNA methyltransferase levels, TMB, and MSI in multiple cancers. GSEA 
also revealed that METTLI6 is involved in the regulation of tumor and immune-related signaling 
pathways. 

[Conclusion] METTL16 is closely related to tumor prognosis and tumor immunoinfiltration, and 
plays an important role in the occurrence and development of various tumors. METTL16 may be a 
promising prognostic biomarker for various tumors and provide new ideas for future tumor 
immune research and therapeutic strategies. 
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N6-HFH ÆRE CN6-methyladenosine, m6A) 修饰 是 最 常见 和 最 丰富 的 RNA 修饰 方式 ， 

尤其 是 在 真 核 信 使 RNA MRNA) HI, 已 有 研究 表明 m6A 修饰 在 调节 RNA 加 工 、 剪 接 、 
成 核 、 翻 译 和 稳定 性 方面 起 着 至 关 重 要 的 作用 , 并 参与 包括 癌症 在 内 的 多 种 疾病 的 发 病 机 秆 
Ul, 其中, 甲 基 转 移 酶 样 蛋 白 16(Methyltransferase-like16,METTL16) 是 新 发 现 的 m6A 甲 基 转 
移 酶 , 可 催化 某 些 编码 和 非 编码 RNA 发 生 甲 基 化 修饰 中, METTL16 的 异常 表达 和 多 种 癌症 
的 产生 及 预后 有 一 定 的 相关 性 ， 例 如 肝癌 四 、 胃 癌 外 、 乳 腺 癌 钻 、 结 肠 癌 但 、 肺 癌 中 等 , (HE 
前 尚 无 研究 对 泛 癌 中 METTL16 的 作用 机 制 进行 全 面 研究 。 
丸 此 ， 本 研究 通过 从 TCGA 数据 库 及 TIMER 等 数据 库 获得 METTLI6 表达 与 泛 癌 之 间 
的 相关 数据 进行 可 视 化 (Kruskal-Wallis 检验 分 析 、Kaplan-Meier 生存 分 析 、 单 因素 Cox [nl 
归 分 析 、GSEA 分 析 )， 分 析 METTL16 在 各 个 肿瘤 中 的 表达 ， 以 及 其 与 肿瘤 预后 、 和 免疫 微 
环境 、 基 因 突 变 和 信号 通路 之 间 的 关联 ， 挖 据 METTL16 在 泛 癌 中 的 潜在 价值 并 加 深 对 它 影 
响 机 制 的 了 解 。 
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1 方法 
1.1 数据 采集 


首先 ， 根 据 实验 要 求 ， 我 们 选取 TCGA 数 据 Æ (The Cancer Genome Atlas, TCGA) 
(https: //portal.gdc.cancer.gov/) 和 基因 型 与 基因 表达 量 关 联 数据 库 (Genotype-Tissue 
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Expression, GTEx (https: //gtexportal.org/) 两 个 数据 库 作 为 了 解 各 个 肿瘤 及 正常 组 织 中 的 
基因 集 的 数据 来 源 。 其 中 ， 我 们 从 GTEx 数据 集 分 析 中 观察 到 METTL16 基因 在 31 个 组 
织 的 表达 水 平 ， 从 TCGA 数据 库 中 获取 了 基因 在 各 个 肿瘤 样本 中 癌 与 癌 旁 的 差异 ， 整 合 


GTEx 数据 库 中 的 正常 组 织 的 数据 和 TCGA 肿 ; 


座 组 织 的 数据 分 析 27 种 肿瘤 的 表达 差异 。 


然后 ， 使 用 Kruskal-Wallis $5354 9r, METTLI6 在 不 同 正常 组 织 和 不 同 肿瘤 组 织 中 的 表达 


水 平 。* 表示 p<0.05,** 表示 p<0.01，*** 表示 p<0.001。 


1.2METTL16 在 各 种 肿瘤 中 的 预后 分 析 
的 表达 谱 数据 的 基础 上 ,我 们 选取 Kaplan-Meier 生存 分 析 、 单 因素 Cox 


在 TCGA 的 基因 


回归 分 析 两 种 
种 肿瘤 


的 表达 不 
1.3 METTL16 在 各 肿瘤 


方法 ， 


的 免疫 温润 分 析 


研究 METTL16 的 表达 对 肿瘤 预后 的 影响 以 及 METTL16 在 TCGA 33 
上 生存 患者 之 间 的 相关 性 。 


我 们 从 肿瘤 免疫 数据 库 (Tumor IMmune Estimation Resource, TIMER) 下 载 了 33 种 癌症 的 
六 种 免疫 浸润 细胞 的 得 分 数据 , 然后 分 别 分 析 METTL16 基因 表达 与 这 些 癌症 免疫 细胞 得 分 
次 ,使 用 R 软件 包 estimate 分 析 了 各 个 肿瘤 样本 的 免疫 评分 和 基质 评分 ， 观 察 


的 相关 性 。 其 


在 33 个 肿瘤 ! 


METTL16 表达 与 免疫 评分 的 相关 性 。 


1.4 METTL16 与 各 肿瘤 免疫 检查 点 之 间 的 相关 性 
我 们 从 UCSC(https://xenabrowser.net/) 数 据 库 中 下 载 了 各 个 肿瘤 的 表达 数据 集 ,分 析 


METTL16 基因 及 四 十 多 个 免疫 检查 点 基 


五 类 免疫 途径 


我 们 使 


的 标 ? 


关 性 是 显著 的 而 且 是 正 相 关 。 
1.6 METTL16 5j TMB 和 MSI 的 关系 


肿瘤 突变 负荷 (tumor mutational burden; TMB), JHK Jx Bi rtg Z 


目 ， 是 一 种 可 定量 生物 标志 物 ， 研 究 证 实 TMB 与 免疫 治疗 预后 具有 机 


广泛 存在 于 原核 和 真 核 生 物 基 因 组 


因 在 各 个 样本 中 的 表达 数据 ,并 计算 了 METTL16 和 
己基 因 的 相关 性 。* 表示 显著 相关 p«0.05, ** 表示 显著 相关 p«0.01, *** 

表示 显著 性 格 p<0.001。 
1.5 METTL16 与 MMR 基因 突变 水 平和 DNA 甲 基 转 移 酶 的 相关 性 
] TCGA 的 表达 谱 数 据 ， 研 究 METTL16 和 5 个 MMRS 基因 
并 构建 METTL16 表达 和 四 个 甲 基 转 移 酶 的 可 视 化 图 谱 。 当 P«0.05 和 R>0.20 时 ， 认 为 相 


的 突变 的 相关 性 ， 


胞 中 所 含有 的 突变 数 
目 关 性 。 


微 卫星 DNA 


1， 微 卫星 不 稳定 性 (microsatellite instability,MSD) 可 作为 


评估 肿瘤 预后 的 一 个 指标 。 在 这 里 我 们 分 别 统 计 了 每 个 肿瘤 样本 的 TMB， 分 析 基 因 表达 与 


TMB 的 关系 。 除 此 之 外 ， 我 们 还 分 析 了 基因 表达 与 MSI 的 相关 性 。 其 中 我 们 使 


曼 等 级 相关 系 


B 
o 


1.7METTL16 5j KEGG fll HALLMARK 通路 的 关系 
基因 表达 将 样本 分 成 高 低 两 组 , 使 用 GSEA 分 析 高 表达 组 和 低 表 达 组 的 KEGG 


我 们 根据 


M HALLMARK 通 } 


H, 


H 


LER 


RIER, RU GSEA， 得 出 最 显著 的 前 三 个 通路 可 视 化。 
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2 结果 


2.1 METTL16 在 不 同 肿瘤 中 的 表达 


A 


Kaniak "Wallis test p=0 


B 


W n wy , m) ! inh 


图 1 METTLI6 在 组 织 中 的 表达 水 平 


Figure 1 The expression level of METTL16 in tissues 


注 : A 为 METTL16 
中 癌 组 织 与 癌 旁 组 织 


在 GTEx 数据 库 中 的 表达 水 平 。B 为 METTL16 在 TCGA 数据 库 
! 的 表达 水 平 。C 为 整合 GTEx 数据 库 与 TCGA 数据 库 中 的 数据 后 ， 


METTL16 的 表达 水 平 与 正常 组 相 比 。* 表示 p<0.05,** 表示 p<0.01，*** 表示 p<0.001。 


首先 ， 我 们 应 用 


GTEx 数据 库 评 估 31 个 正常 组 织 中 的 METLL16 表达 。 结 果 表 明 ， 


在 所 有 评估 的 正常 组 织 中 ，METLL16 在 子宫 颈 ， "T3 垂体 ， 子 襄 的 表达 水 平 最 高 ， 在 血 
液 , 胰腺 ， 肝脏 的 表达 水 平 最 低 。(Kruskal-Wallis 检验 , p=0) (图 1A)。 此 外 , 基于 TCGA 
数据 库 ， 我 们 发 现 METLL16 在 BLCA, BRCA, KICH, LUAD, LUSC, ULEC, THCA, 
PRAD 的 表达 水 平 相 对 较 高 (图 1B), Æ CHOL, COAD, HNSC, LIHC, STAD, ESCA, 


KIRP,LGG 的 表达 水 


FE 相对 较 低 (图 1B)。 考 虑 到 TCGA 中 缺乏 正常 样本 对 照 的 肿瘤 ， 我 


们 将 来 自 GTEx 数据 库 的 相应 正常 样本 与 TCGA 数据 库 整 合 为 图 1C， 来 评估 不 同 癌 症 中 


METTL16 表达 的 差异 。 结 果 表 明 ，METLL16 在 CHOL、ESCA、GBM、HNSC、KIRP、 


LAML, LGG, PAAD, 


LUAD, LUSC. OV. 


SKCM, STAD, TGCT 表达 上 调 (图 1C), 在 BLCA、CESC、COAD、 
PRAD, THCA. UCEC, UCS 表达 下 调 (图 1C). 
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2.2 METTL16 在 不 同 肿瘤 中 的 预后 分 析 


Bm f Á— BLCA — 1000994100) 
i «i| 1 BRCA — 1030014105) 
TT CESC — 0909-109) 

H OPETE NES. | i marmon CHOL 105095 1.17) 
142-2 2 —1— ái 1cz-2—2-2 12] " "am 
i ; i T DUC — 0908-110) 
c ESA — 081082100) 
GB) 1010985104 

HNSC — 099096102) 

| KICH 120.01 141) 
i KRC — 094091-097) 
KAP — 099094100) 

n UM, — 099095104) 
¥ M6 — 1009-19) 
E URC 104098 1.11) 
E WAD 1009710) 
USC 102097104) 

MES — 099094104) 

i O 4010985103 
H PAD  099088=098) 
RPG 097081118) 

s. PRAD — 104082110) 
i RED 0890791) 
H SC — 10010) 
a SCU — 1009-10) 
SWO — 09908510) 
TOT MNMS1U) 

j| THCA 101(094*1.1) 
i TM — 0807-09) 
WEG — (9610) 

$ WS — 108096710 
t UW — 1090005109) 


P Value 
DD. eses 3$--—- 


0n 10 1a 
HR(SSCI 


2 METTLI6 表达 与 癌症 患者 的 总 体 生 存 率 的 关系 


Figure 2 Relationship between METTL16 expression and overall survival of cancer patients 


注 : 


A 为 METTL16 高 低 水 平 与 BRCA OS 之 间 关 系 的 K-M 生存 曲线 ， 
B 为 METTL16 高 低 水 平 与 KICH OS 之 间 关 系 的 K-M 生存 曲线 ， 
C 为 METTLI6 高 低 水 平 与 KIRC OS 之 间 关 系 的 K-M 生存 曲线 ， 
D 为 METTL16 高 低 水 平 与 PAAD OS 之 间 关 系 的 K-M 生存 曲线 ， 
E 为 METTL16 高 低 水 平 与 READ OS 之 间 关 系 的 K-M 生存 曲线 ， 
F 为 METTL16 高 低 水 平 与 SARC OS 之 间 关 系 的 K-M 生存 曲线 ， 
G 为 METTL16 高 低 水 平 与 THYM OS 之 间 关 系 的 K-M 生存 曲线 ， 
H 为 METTL16 高 低 水 平 与 UAM OS 之 间 关 系 的 K-M 生存 曲线 ， 


(BRCA: 乳 腺 浸润 癌 :KICH: 肾 嫌 色 细胞 癌 ; KIRC: 


肾 透 明细 胞 癌 ; PAAD: 胰腺 癌 ; READ: 


直肠 腺 癌 ; SARC: 肉瘤 ; THYM: 胸腺 癌 ; UVM: 葡萄 膜 黑 色素 瘤 )。 


METTL16 与 肿瘤 预后 的 关系 分 析 结 果 显 示 ， 
KICH (HR=1.2, P=0.03300)、 KIRC (HR=0.94, P=0 


UVM (HR=1.05，P=0.00560〉 中 存在 预后 意义 (图 


METTL16 在 BRCA(HR=1.03, P=0.014)、 


.00021)、PAAD (HR=0.93, P=0.00520)、 
READ (HR=0.89, P=0.0440)、 SARC CHR-1.02, P-0. 


03700), THYM (HR=0.88, P=0.04100)、 
2D. K-M 分 析 曲 线 结 果 显 示 ，BRCA、 


KICH.SARC.UVM 中 METTL16 高 表达 水 平 以 及 KIRC、PAAD、READ THYM 中 METTL16 


低 表 达 水 平 与 不 良 预 后 有 关 ( 图 2 A-H)。 
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2.3 不 同 肿瘤 中 METTL16 与 肿瘤 免疫 浸润 和 肿瘤 微 环境 之 间 的 相关 性 


ioo2(METTL1G Tee 1) 


Timer Database score. 


ESTIMATESGre 


oni 


图 3 METTL16 表达 与 不 同类 型 癌症 中 的 免疫 浸润 水 平 的 关系 


EEL] 
emt 


Figure 3 Relationship between METTL16 expression and the level of immunoinvasion in different 


types of cancer 


iX: A 为 METTL16 的 表达 水 平 与 COAD、KIRC、LIHC 中 免疫 细胞 浸润 之 间 的 相关 性 分 


析 ，B 为 METTL16 的 表达 水 平 与 免疫 评分 ， 基 质 评分 以 及 ESTIMATE 的 免疫 订 


系 o 


表 1 METTLI6 KASR JE 


分 的 关 


分 以 及 ESTIMATE 的 免疫 评分 关系 


Table 1 The relationship between METTL16 expression and immune score, matrix score and 
immune score of ESTIMATE 


项 目 


免疫 评分 SARC 


iz 


E VEA SARC 


ESTIMATE 免疫 评分 SARC 


COAD 


METTL16 
I 值 P 值 
0.363 <0.001 
0.267 <0.001 
0.355 <0.001 
-0.363 <0.001 
-0.192 <0.001 
-0.047 <0.001 
-0.363 <0.001 
-0.192 <0.001 
0.267 <0.001 


ik: SARC, AX; COAD， 结 肠 癌 ; PAAD， 上 胰腺 ; THCA, MORIRE; BRCA， 乳 腺 浸润 


JA 
Th o 
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通过 对 METTL16 与 免疫 细胞 浸润 以 及 肿瘤 微 环 境 之 间 的 相关 性 进行 评分 ， 有 助 于 预 
测 癌症 名 的 预后 和 治疗 敏感 性 。 首 先 ， 我 们 利用 TIMER 数据 库 来 研究 METTL16 与 人 类 
33 种 癌症 的 六 种 免疫 温润 细胞 之 间 的 潜在 关联 ,结果 显示 , METTL16 的 表达 水 平 与 COAD、 

KIRC、LIHC 等 肿瘤 的 免疫 浸润 水 平 的 相关 性 最 强 。 而 且 METTL16 的 表达 水 平均 与 COAD、 
KIRC、LIHC 中 的 六 种 免疫 浸润 细胞 即 B cH. CD4- T 细胞 、CD8+TI 细胞 、 树 突 状 细 
胞 、 巨 哈 细 胞 和 中 性 粒 细 胞 的 水 平 呈正 相关 (图 3A)。 此 外 ,我 们 整合 了 免疫 评分 、 基 质 评 
分 和 ESTIMATE 免疫 评分 , 以 进一步 评估 METTL16 表达 与 癌症 免疫 浸润 之 间 的 关系 。 结 
果 表 明 ，METTL16 的 表达 与 SARC (R-0.363,P«0.001 )、COAD(R=0.267,P<0.001)、 
PAAD(R-0.355,P«0.001) 的 免疫 评分 均 呈 正 相 关 。 METTLI6 的 表达 与 
SARC(R--0.363,P«0.001). THCA(R--0.192,P«0.001). BRCA(R--0.047,P«0.001) f] & Ji VE 4] 
均 呈 负 相 关 。METTL16 的 表达 与 COAD(R-0.267,P«0.001) ESTIMATE 的 免疫 评分 呈正 相 
X. 5 SARC(R--0.363,P«0.001) ~ THCA(R--0.192,P«0.001) 呈 负 相关 (图 3B、 表 1). 


2.5 不 同 肿瘤 中 METTL16 的 表达 水 平 与 MMR 基因 突变 水 平和 DNA 甲 基 转 移 酶 的 关 
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5 METTL16 在 泛 癌 中 的 表达 与 5 个 MMR 基因 和 4 个 DNA 甲 基 转移 酶 的 相关 性 
Figure 5 The correlations of METTL16 expression with five MMR genes and four DNA 
methyltransferases in pan-carcinomas. 
iE: CA) METTL16 的 表达 水 平 与 癌症 中 五 种 MMR 基因 表达 水 平 的 相关 性 ， 
(B) METTL16 的 表达 水 平 与 癌症 中 四 种 DNA. 甲 基 转移 酶 的 相关 性 ， 
* 表 示 相 关 性 差异 为 P<0.01，*** 表 示 相 关 性 差异 为 P<0.001。 


为 了 确定 METTL16 在 肿瘤 进展 中 的 作用 ， 我 们 评估 了 METTL16 表达 水 平 与 MLH1、 
MSH2、MSH6、PMS2、EPCAM 五 个 MMR 基因 突变 水 平 的 关联 ， 结 果 显 示 ， 除 了 DLBC 
之 外 ， 各 个 肿瘤 中 METTL16 的 表达 水 平 与 5 个 MMR 基因 突变 水 平 显著 相关 ， 且 几乎 均 
为 正 相 关 。 (图 SA)。 除 此 之 外 ， 我 们 评估 了 METTL16 表达 水 平 与 DNMT1、DNMT2、 
DNMT3A 、DNMT3B 四 种 DNA. 甲 基 转 移 酶 之 间 的 关系 ,结果 显示 , 除了 DLBC 和 UCS 
之 外 ， 各 个 肿瘤 中 METTL16 的 表达 水 平 与 四 种 DNA 甲 基 转 移 酶 的 水 平 几乎 均 呈 显著 正 
相关 CP «0.05,R >0.20) (图 5B). 
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2.6 不 同 肿瘤 METTL16 表达 水 平 与 TMB 和 MSI 的 关系 
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7 METTL16 表达 与 肿瘤 突变 负荷 CTMB )、 微 卫星 不 稳定 性 (MSI 的 关系 
Fig7 Correlation between METTL16 expression and tumor mutation burden (TMB), 
microsatellite instability (MSI). 


METTLI6 表达 与 TMB 的 关系 图 显示 ，METTL16 在 COAD (了 P<0.0001) ~ 
LAML(P=0.032)、 LGG(P=0.00032)、MESO(P<0.0001)、UCEC(P<0.0001) 中 的 表达 与 TMB 
呈正 相关 ， 在 BLCA(P=0.032)、BRCA(P<0.0001)、CESC(P=0.032)、ESCA(P=0.00023)、 
KIRC(P-0.039), PRAD(P-0.024), PRAD(P-0.028), THCA(P«0.0001), UVM(P-0.033) 
与 TMB 呈 负 相关 (图 7A)。 此 外 ，METTL16 表达 与 MSI 的 关系 图 显示 ，METTL16 在 
CHOL(P=0.43)、COAD(P<0.0001)、LUSC(P<0.0001)、UCEC(P=0.0086) 中 的 表达 与 MSI 
呈正 相关 ， 在 PRAD(P<0.0001)、SKCM(P=0.0071)、THCA(P=0.00017 与 MSI 呈 负 相关 
(图 7B). 


2.7METTL16 在 KEGG 和 HALLMARK 通路 的 GSEA 分 析 
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图 6 METTLI6 的 KEGG fll HALLMARK 通路 富 集 图 
Fig 6 KEGG and HALLMARK pathway enrichment of METTL16 


我 们 通过 GSEA 研究 METTL16 的 生物 学 作用 。KEGG 分 析 显 示 METTLI16 在 长 时 程 
增强 (P=0)、 赖 氨 酸 退化 (P=0) 和 WNT 信号 通路 (P=0) 中 高 表达 。 此 外 ，HALLMARK 分 
析 显 示 METTL16 在 血红 素 代谢 (P=0.002)、 响 应 紫外 线 辐射 (P=0) 和 和 蛋白质 分 泌 (P=0) 
通路 中 高 表达 (图 6) 。 


3 讨论 
泛 癌 分 析 是 指 利用 多 个 数据 库 分 析 不 同 肿瘤 中 基因 的 表达 、 预 后 、 突 变 模式 和 功能 的 生 
物 信 息 学 分 析 ， 它 可 以 揭示 不 同 肿 瘤 之 间 的 异同 ,对 癌症 的 预防 , 诊断 和 治疗 以 及 指导 癌症 
相关 研究 具有 重要 意义 匀 。N6- 腺 苷 甲 基 化 〈N6 - methyladenosine, m6A ) 是 一 种 高 度 保守 的 
内 部 RNA 修饰 ， 可 以 影响 影响 RNA 加 工 的 多 个 方面 ， 包 括 空间 结构 的 折合、 稳定 性 、 半 
衰 期 控制 、 前 体 mRNA 的 剪接 、 输 出 和 翻译 0。METTL16 是 m6A RNA 甲 基 转 移 酶 中 的 
一 种 ， 已 有 多 项 研究 表明 其 与 胃癌 00、 肺 癌 四 、 前 列 腺 癌 59、 卵 梨 癌 03] 等 多 种 癌症 的 发 生 
| 及 其 预后 相关 。 为 进一步 研究 METTL16 与 各 种 癌症 的 关系 , 本 研究 整合 了 TCGA 与 GTEx 
~ 数据 库 的 样本 从 基因 表达 、 免 疫 分 析 等 角度 对 其 进行 了 一 系列 的 泛 癌 分 析 。 
. 我 们 的 研究 显示 ,METTL16 对 不 同 癌症 的 预后 有 不 同 影响 ,同时 我 们 发 现 METTL16 的 
表达 水 平 与 结肠 癌 (COAD)、 肾 透明 细胞 癌 (KIRC)、 肝 癌 (IHC) 中 的 的 六 种 免疫 浸 
润 细胞 即 B S. CD4- T 细胞 、CD8+T 细胞 、 树 突 状 细胞 、 巨 唉 细胞 和 中 性 粒 细胞 的 
e 水 平 呈正 相关 。 此 外 ,METTL16 表达 也 与 多 种 癌症 中 的 TMB/MSIDNMTS/MMR 基因 相关 。 
T GSEA 分 析 的 结果 显示 METTL16 在 长 时 程 增强 、 赖 氨 酸 退化 、WNT 信号 通路 、 血 红 素 代 
本 谢 、 响 应 紫外 线 辐射 和 蛋白 质 分 泌 通 路 中 高 表达 。 
Dai YZ 等 人 的 实验 结果 表明 METTL16 在 肝癌 CLICH/HCCO 中 的 表达 水 平 上 调 ， 与 
我 们 的 研究 结果 一 致 叫 。 此 外 我 们 的 研究 还 发 现 METLLI6 在 胆管 癌 (CHOL)、 结 肠 癌 
CCOAD )、 头 颈 癌 CHNSCO. 胃癌 (STAD)、 食 管 癌 (ESCA)、 肾 乳头 状 细 胞 癌 (KIRP)、 
脑 低级 别 胶 质 瘤 〈LGG ) 中 的 表达 水 平 相对 较 高 ， 而 在 膀胱 尿 路 上 皮 癌 〈BLCA)、 乳 腺 浸 
润 癌 (BRCA )， 肾 嫌 色 细胞 癌 (KICH)、 肺 腺 癌 (LUAD)， 肺 鳞 状 细胞 癌 (LUSC)， 子 宫 内 膜 
本 羊 癌 (UCEC), 甲状 腺 癌 (THCA), 胰腺 癌 (PAAD ) 的 表达 水 平 相对 较 低 。 表 明 , 今后 METTL16 
- 可 成 为 上 述 癌症 中 预后 观察 的 指标 ， 但 仍 需 进一步 实验 研究 验证 。 
FE Wang XK 等 人 的 研究 表明 METTLI16 在 胃癌 具有 促 癌 作用 ， 同 时 METTLI6 的 表达 水 
平 与 胃癌 的 不 良 预后 呈正 相关 093。 但 是 我 们 的 研究 结果 表明 METTL16 对 胃癌 的 预后 影响 较 
小 ， 可 能 需要 更 大 样本 量 实验 研究 进行 验证 。 此 外 ， 我 们 研究 结果 显示 ，METTL16 高 表达 
水 平 与 乳腺 浸润 癌 (BRCA)、 肾 嫌 色 细胞 癌 (KICH)、 肉 瘤 (SARC)、 眼 部 黑色 素 瘤 CUVM) 
不 良 的 预后 有 关 ， 而 与 肾 嫌 色 细 胞 癌 (KICH), ARE (PAAD, EU (READ), Hie 
(THYM) 较 好 的 预后 有 关 。 
WNT 信号 通路 具有 高 度 的 保守 型 ， 研 究 表 明 其 对 胚胎 发 育 、 细 胞 周期 调节 、 炎 症 和 癌 
症 在 内 的 许多 过 程 有 重要 影响 , 当 WNT 及 其 下 游 效应 器 在 失控 时 会 促进 肿瘤 的 发 生 、 生 长 
和 转移 09。 既 往 研 究 显示 ， 多 种 癌症 中 存在 WNT 信和 号 通路 的 异常 表达 ， 例 如 肝癌 07I9、 结 
直肠 癌 、 皮 肤 癌 09、 乳 腺 癌 P,20 等 。 
WNT 信号 通路 分 为 经 典 (B-catenin 依赖 性 ) 和 非 经 典 (B-catenin 非 依赖 性 ) 两 种 ,其 中 以 前 
者 的 研究 较为 广泛 。 研 究 表 明 ， 经 典 (B-catenin 依赖 性 JWNT 通路 的 基因 突变 或 受 体 激活 与 
组 织 的 癌变 有 很 大 的 关联 性 , 其 可 通过 控制 细胞 的 异常 生长 和 凋 亡 的 进程 、 诱 导 癌 细胞 转移 
和 粘 附 、 加 速 肿瘤 血管 的 生成 , 促进 癌 细 胞 的 自我 更 新 等 机 制 影响 癌症 的 发 生 和 预后 , 此 外 ， 
还 有 研究 表明 经 典 WNT 通 路 子 在 免疫 工 细胞 的 发 育成 熟 以 及 癌症 免疫 逃逸 发 挥 重要 的 作用 
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04。 然 而 ， 非 经 典 的 WNT 通路， 最 近 研 究 表 明 ， 其 Wnt/Ca2+ 途 径 与 造血 干细胞 的 衰老 、 
破 骨 细胞 的 分 化 有 关 ， 同 时 也 影响 着 某 些 癌症 的 发 展 扬 约 。KEGG 的 分 析 显 示 METTLI6 在 
WNT 信号 通路 (P=0) 中 高 表达 。Xiao-Kun Wang 等 人 的 研究 发 现 ，METTL16 可 以 通过 上 
调 胃 癌 细 胞 周期 蛋白 D1 的 表达 来 促进 胃癌 细胞 的 增殖 03。DAVID TAI 等 人 的 研究 表明 , 在 
B -连环 蛋白 依赖 的 途径 His WNT 刺激 的 缺失 会 导致 8 -连环 蛋白 的 破坏 ， 致 使 细胞 周期 蛋 
EDI 表达 水 平 上 调 09。 这 些 研 究 表明 METTL16 对 癌症 的 发 生 及 其 预后 的 影响 可 能 与 WN 
通路 的 高 表达 状态 有 一 定 的 关系 ， 然而 具体 的 机 制 还 需要 进一步 的 研究 。 

本 课题 还 针对 METTL16 的 表达 水 平和 TMB、MSI、MMR、DNTSR 的 关系 进行 了 相 
的 研究 。 肿 瘤 突变 负荷 (MB) 和 微 卫 星 不 稳定 性 MSD 被 认为 是 影响 肿瘤 发 生 和 进展 的 
重要 因素 。 肿 瘤 突变 负荷 (tumour mutation burden, TMB) 是 一 项 肿瘤 预测 指标 ， 其 与 基因 
突变 有 关 ， 包 括 点 突变 和 插入 缺失 鸣 。 同 时 ， 微 卫星 不 稳定 (microsatellite instability, MSD 也 
和 DNA HX, Æ DNA 复制 过 程 中 ， 串 联 重 复 区 域 易 发 生 DNA 复制 错误 P5。 因 错 配 修复 
(mismatch repair, MMR) 蛋白 可 以 起 到 修复 DNA 的 复制 错误 的 序列 的 作用 ，MMR 基因 的 突 
变 或 甲 基 化 均 可 引起 MSI 的 发 生 P9。 本 研究 结果 显示 ,METTL16 表达 水 平 与 MLH1、MSH2、 
MSH6、PMS2、EPCAM 五 个 MMR 基因 突变 水 平 的 关联 及 四 种 DNA. 甲 基 转 移 酶 的 水 平 
均 呈 正 相 关 , 且 具有 一 定 统计 学 意义 。 此 外 ,METTL16 表达 与 TMB 的 关系 图 显示 , METTLI6 
在 COAD、UCEC 中 的 表达 与 TMB 呈正 相关 ， 在 PRAD、THCA 5 TMB 呈 负 相关 。 此 外 ， 
METTL16 表达 与 MSI 的 关系 图 显示 ，METTL16 在 COAD、UCEC 中 的 表达 与 MSI 呈正 相 
关 ， 在 PRAD、THCA 与 MSI 呈 负 相关 。 

虽然 本 次 研究 通过 生 信 分 析 确 定 了 METTL16 在 泛 癌 中 的 潜在 价值 ， 但 由 于 本 次 研究 
只 是 基于 数据 库 的 分 析 得 出 的 结果 ,没有 经 过 临床 样本 测试 进行 验证 ， 故 还 存在 一 定 的 局 限 
性 。 未 来 的 研究 可 以 朝 此 方向 努力 ， 同 时 也 可 以 更 加 深入 的 发 气 METTL16 在 肿瘤 发 生 和 
发 展 过 程 中 的 具体 机 制 ， 丰 富 METTL16 基因 的 临床 价值 。 
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